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Nonstandard Analysis and Shock Wave Jump Conditions in 
a One-Dimensional Compressible Gas 

by 

Roy S. Baty 
F. Farassat 
and 

John A. Hargreaves 

ABSTRACT 

Nonstandard analysis is a relatively new area of mathematics in which infinitesimal 
numbers can be defined and manipulated rigorously like real numbers. This report presents a 
fairly comprehensive tutorial on nonstandard analysis for physicists and engineers with many 
examples applicable to generalized functions. To demonstrate the power of the subject, the 
problem of shock wave jump conditions is studied for a one-dimensional compressible gas. ft is 
assumed that the shock thickness occurs on an infinitesimal interval and the jump functions in the 
thermodynamic and fluid dynamic parameters occur smoothly across this interval. To use 
conservations laws, smooth pre-distributions of the Dirac delta measure are applied whose 
supports are contained within the shock thickness. Furthermore, smooth pre-distributions of the 
Heaviside function are applied which vary from zero to one across the shock wave. It is shown 
that if the equations of motion are expressed in nonconservative form then the relationships 
between the jump functions for the flow parameters may be found unambiguously. The analysis 
yields the classical Rankine-Hugoniot jump conditions for an inviscid shock wave. Moreover, 
non-monotonic entropy jump conditions are obtained for both inviscid and viscous flows. The 
report shows that products of generalized functions may be defined consistently using 
nonstandard analysis; however, physically meaningful products of generalized functions must be 
determined from the physics of the problem and not the mathematical form of the governing 
equations. 
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1.0 INTRODUCTION 


The lion’s share of results in modem mathematics have been discovered from the study 
and application of partial differential equations. Currently, no known mathematical theory 
provides techniques to study all the linear and nonlinear partial differential equations of 
mathematical physics. The lack of such a coherent theory makes the analysis of abstract linear 
partial differential equations very difficult, and the analysis of general nonlinear partial 
differential equations essentially insurmountable. The effect of the absence of a central theory is 
that each different type of physical problem, and hence, each type of equation, requires a different 
analytical or numerical method for its solution. Therefore, the subject of partial differential 
equations reduces to a potpourri of seemingly unrelated techniques. 

Existing methods to analyze partial differential equations are not very useful for nonlinear 
problems, because almost all methods used are based on the linear mathematical structures 
developed in functional analysis. Since linear mathematical structures do not describe the 
simplest nonlinearity (multiplication) it is not surprising that functional analysis has not yielded a 
general framework to study partial differential equations. 

Many approaches have been used to attempt to define a useful multiplication operation on 
spaces of functionals applied in linear analysis to study differential equations. The motivation for 
multiplying functionals comes from a plethora of applications that occur in modeling physical 
problems: construction of solutions of certain partial differential equations requires the 
multiplication of linear functionals. However, a fundamental theorem of Schwartz [1] shows that 
the standard algebraic product cannot be defined consistently on the spaces of functionals, such as 
the distributions or classical generalized functions; all such methods will produce ambiguous 
results leading to logical and arithmetical contradictions. Since standard multiplication does not 
produce meaningful results, products of generalized functions have been defined which are 
weaker than standard multiplication. The resulting spaces of generalized functions are algebras 
containing the classical generalized functions as a proper subset. 

Over the last thirty years, E.E. Rosinger [2-5] and J.F. Colombeau [6,7] working 
independently, have developed a theoretical framework for quotient algebras of generalized 
functions to analyze nonlinear partial differential equations. Their approach extends classical 
linear mathematics to nonlinear mathematics by defining multiplication on the classical 
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generalized functions by embedding these linear spaces into abstract algebras of nonlinear 
generalized functions with a weak operation of multiplication. The potential use of such theories 
of nonlinear generalized functions allowing the analysis and approximation of solutions to 
nonlinear partial differential equations is enormous. Such a theory could impact almost all areas 
of science and technology. In engineering, for example, a coherent mathematical theory would 
help advance the numerical methods used to predict solutions of complicated problems and could 
reduce significantly the physical experimentation needed in the solution of engineering problems. 

While it is probable that these theories of generalized functions will yield insight into the 
theory and computation of certain classes of nonlinear field equations, significant potential 
limitations of the theory of nonlinear generalized functions are known. Since any logically 
consistent operation of multiplication of generalized functions must be much weaker than the 
classical multiplication of numbers and functions, such a product of generalized functions will 
always yield large collections of abstract objects representing weak solutions of the problem 
under investigation. The upshot is that a theory of nonlinear generalized functions allows the 
solution of problems not defined in classical mathematics by relaxing the uniqueness of 
multiplication. Basically, problems with no classical weak solutions, such as the product 8 • H , 
are converted into problems with uncountably many generalized solutions. The central issue then 
becomes how to determine a physically observed solution from an infinitude of possible weak 
solutions. 

An example of both the power and limitation of the theory of nonlinear generalized 
functions can be seen in the research of Rosinger [5,8], In his work, conditions are given which 
yield generalized solutions for arbitrary polynomial nonlinear partial differential equations. 
Rosinger has shown that generalized solutions can be obtained for all equations of the form: 


a L|3(a) 


= /O) 


( 1 . 1 ) 


where U is an unknown function d is a derivative operator, a, p are multi-indices, and c,f are 
specified functions. However, this result is not immediately useful in practice, because it tells 
nothing about how to approximate the solution U . 

To approximate solutions, Colombeau (see references [6,7]) constructed an explicit theory 
of nonlinear generalized functions. This theory has been applied to analyze and compute solutions 
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to nonlinear partial differential equations. The logic of solving a differential equation in a space of 
Colombeau generalized functions follows the standard approach of functional analysis: a method 
is established for generating sequences to approximate solutions within a collection of nonlinear 
generalized functions and then a compactness argument is applied to show that the approximating 
sequences actually converge to limits within the space of approximating objects. A resulting limit 
then represents a solution to the equation. A key difficulty in applying a space of Colombeau 
generalized functions is determining a specific representation for the products of generalized 
functions occurring in the differential equation. As Colombeau [9] and his co-workers [10,11] 
have shown, different physical problems require different distinct representations for the same 
product of generalized functions such as 8 • H . The representation of the product is found by 
adding more physical information about the problem. Roughly speaking, the solution of a 
concrete problem via nonlinear generalized functions requires: the equations governing the 
physics, the initial-boundary values, and additional problem-specific physical information 
required to fix the product of generalized functions. For example, in continuum shock wave 
problems, the additional information is obtained by requiring the product of generalized functions 
to satisfy the standard shock wave jump conditions. 

In studying the Colombeau generalized functions a number of authors including, Todorov 
[12], Oberguggenberger [13], and Hoskins and Sousa Pinto [14], have shown that versions of the 
Colombeau algebras of nonlinear generalized functions may be constructed using nonstandard 
analysis. Nonstandard analysis is a relatively new area of mathematics which emerged from basic 
research in mathematical logic. The subject was discovered in the early 1960s by A. Robinson 
[15]. The main contribution of nonstandard analysis to mathematics is the extension of the real 
numbers R to the hyperreal numbers *R which contain infinitely small (infinitesimal) and 
infinitely large numbers. The infinitesimals in the hyperreal number system have the properties 
that Leibniz applied in his development of the calculus, which justifies the algebraic 
manipulations of infinitesimals that engineers and physicists often use. Moreover, the existence of 
distinct infinitely large numbers that can be manipulated and used like finite numbers to solve 
problems provides a powerful new tool for applications. Nonstandard methods may be used to 
construct algebras of nonlinear generalized functions or nonstandard analysis may be used 
directly to study problems involving multiplication of generalized functions. The second 
approach is considered here. 
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In this report the problem of shock wave jump conditions for a one-dimensional 
compressible gas is studied using nonstandard analysis. A class of nonstandard Heaviside 
functions is applied to model the microstructure of the shock wave jump. The term microstructure 
refers to the analytical description of the functions which represent the flow variables such as 
density, pressure, velocity and entropy across the shock wave. The shock wave is governed by the 
standard conservation equations written in primitive form. The shock jump is then modeled to 
occur over an infinitesimal interval. The family of Heaviside functions applied has some very 
interesting properties which are exploited in the subsequent analysis: for a function (j)(x) with 
jump [())] = (j) r - atx = 0 define 

cKx) = ^,+ [») (1.2) 

where cj)(nr) is at least a C 1 function and II(x) is a differentiable, nonstandard Heaviside 
function. Following Colombeau [16], each flow parameter is assumed to have a different 
Heaviside function associated with it at the shock. This means that each flow variable may have a 
different shape or microstructure across the shock wave. 

It is shown that if the various Heaviside functions for the flow parameters have their 
jumps located on the same infinitesimal interval, then the governing equations expressed in 
nonconservative form yield the relations between all these Heaviside functions unambiguously. 
Strictly speaking, the mathematics of the product of nonstandard generalized functions is 
unambiguous; that is, use of differentiable nonstandard Heaviside functions removes the 
mathematical difficulties associated with products such as 8 • H . The micro structure of the 
Heaviside functions predicts the physical microstructure of the shock wave jump if the governing 
equations of motion model the physics of the shock layer. For an inviscid shock wave, this 
analysis cannot reproduce the physical microstructure, because the physical information 
describing the shock region is missing. However, for a viscous shock wave, this analysis 
reproduces the experimentally observed shock structure. The analysis also yields the classical 
Rankine-Hugoniot jump conditions found for inviscid compressible flows. Moreover, the unusual 
microstructure of the Heaviside function for the inviscid entropy jump is obtained by the 
nonstandard method and is similar to that found by Morduchow and Libby [17] using a viscous 
flow description and the results of Iollo and Salas [18] using Colombeau’s nonlinear generalized 
functions. 



The present work was motivated by the authors’ interest (see Baty et al. [19]) in 
understanding how to add physical information to specify nonlinear generalized functions 
contained in differential equations which model problems with well defined regions of solution 
discontinuity, such as a shock wave jump problem. In an intuitive sense Colombeau’s generalized 
functions must represent the infinitesimal microstructure of the shock wave across the jump. 
Therefore in analogy with these generalized functions, the present work applies nonstandard 
analysis to build a rigorous, explicit example of the infinitesimal microstructure of a shock wave. 
If the existence of the nonstandard Heaviside functions is assumed given, the present 
development greatly simplifies the description of the product of generalized functions, because 
the level of mathematical abstraction is reduced to an advanced calculus problem and not that of a 
research problem in functional analysis, Keisler [20] and Laugwitz [21], Moreover, the derivation 
of the hyperreal shock wave jump conditions helps distinguish between mathematical and 
physical ambiguities in the inviscid and viscous flow cases. 

Chapter 2 begins the technical development of the report by summarizing the equations of 
motion in nonconservative form for a compressible perfect gas in one spatial dimension. 
Equations for the inviscid and viscous flow cases as well as the non-heat conducting and heat 
conducting cases are given. The basic physical assumptions about the gas are also outlined. 

Chapter 3 outlines a generic construction of the hyperreal number system. Examples of 
nonstandard differentiable Heaviside functions are then introduced. Some basic properties of the 
nonstandard Heaviside functions are reviewed and examples of products of these nonstandard 
jump functions and their derivatives are exhibited in the context of products of generalized 
functions. 

Chapter 4 applies nonstandard differentiable Heaviside functions to derive shock wave 
jump conditions for an inviscid compressible gas. The relations between the jump functions are 
obtained and then used to show that the entropy must exhibit a peak in the shock layer. The 
location of the entropy peak is found to depend on the Heaviside functions for the specific volume 
and pressure and cannot be specified for an inviscid flow field. 

Chapter 5 extends the derivation of shock wave jump conditions of Chapter 4 to a viscous 
compressible gas and shows that the flow field Heaviside functions are specified by the viscous 
flow assumption. The shock wave jump functions are computed for both the non-heat conducting 
and heat conducting cases for low Mach numbers. The computed Heaviside function for density is 
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also compared with the density measured for a normal shock wave in a low density helium flow 
field. 

Chapter 6 then concludes the report with a summary of the key observations and results 
presented. 
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2.0 EQUATIONS OF MOTION 


2.1 Conservation Equations 

The Navier-Stokes equations govern the motion of compressible fluid flow. In cartesian 

12 3 

coordinates ( x ,x ,x ) , the equations governing conservation of mass, momentum and energy 
are given by: 


|_P + /5p + 

8t dx k 



= o, 


( 2 . 1 ) 



, k 
+ pit 


dii 

dx 


dp_ + _d_ 

dx 7 dx ' 


k\ 
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A. T 

V dx J 


+ . 


dx 
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du + dii 
vdx 7 dx) 


( 2 . 2 ) 


de kde du' d 

8t dx dx dx 


dL 

dx)) 


+ X 


f \ £\ 2 
ou 


Kdx k J 


+ p 


1 du + du ^ 
ISx 7 dx') 


dii 

dx 1 


(2.3) 


k 

where u (k = 1, 2, 3 ) are the velocity components, p is the density, p is the pressure, e is the 
energy per unit mass, and T is the temperature. The parameters p, A and k are the dynamic 
viscosity, the second viscosity coefficient and, the thermal conductivity of the fluid, respectively. 
These parameters are assumed to be constant. The conservation equations (2.1), (2.2) and (2.3) 
represent five equations in terms of seven physical variables. To close the system of equations two 
equations of state must be added which relate the thermodynamic variables p, p, e and T : 


P = P(P, 2), 

(2.4) 

e = e(p, T). 

(2.5) 


Developments of equations (2.1), (2.2), and (2.3) and the equations of state (2.4) and (2.5) may be 
found in Currie [22] and Thompson [23]. 


2.2 Inviscid Equations 

For the development of the shock wave jump conditions for an inviscid, non-heat 
conducting gas the conservation equations are assumed to model the motion of a gas in one spatial 
dimension and time: 


(u , if, u ) = (u, 0, 0) so that u = u and x = x 


(2.6) 
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Furthermore, since the gas is assumed to be inviscid and non-heat conducting the fluid parameters 
drop out: 

p = 0 , X = 0 , and k = 0 . 

With these assumptions, equations (2.1) to (2.3) become 


dp 

8t 


, dp 
+ u— + 
dx 



= 0 , 


(2.7) 



du 

+ p “5 


dp 
dx ’ 


( 2 . 8 ) 


de , de du 

(2.9) 

y dt dx F dx 

For a perfect gas, the equations of state are given by 


> 

II 

(2.10) 

p = p RT, 

(2.11) 


where c v is the specific heat at constant volume (assumed constant) and R is the specific gas 
constant. 

The conservation equations are expressed in terms of specific volume v = - rather than 

P 

the density. Also, the notation for partial differentiation is simplified to the subscript notation, 

dt 

E, x = -^. With these assumptions equations (2.7) - (2.9) become: 


v t - vu x + uv x = 0 , 

(2.12) 

u t + uu x + vp x = 0, 

(2.13) 

e t + ue x + vpu x = 0. 

(2.14) 


Next combining the equations of state with the energy equation (2.14), it may be shown 

that 


p(s t + us x ) = c v (p t + up x + ypu x ) , (2.1 5) 

where 


s = c v log (pv Y ), 


(2.16) 
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is the entropy per unit mass and where y is the ratio of specific heats. Since the gas motion on 
either side of a shock wave is assumed to be isentropic, equation (2.15) reduces to zero along 
streamlines upstream and downstream of a shock wave. However, a shock wave is not isentropic 
process, so across the shock wave equation (2.15) is not zero: 

p(s t + us x ) * 0 , 


which implies 


c v (p t + up x + ypu x ) * 0 . 


Therefore, the equations governing the unsteady, one-dimensional motion of an inviscid, non-heat 
conducting gas across a shock wave are given by: 


v t - vu x + UV X = 0 , 

(2.17) 

u t + uu x + vp x = 0, 

(2.18) 

S = CylOg(pV Y ). 

(2.19) 


2.3 Viscous Equations 

For the development of the shock wave jump conditions for an viscous, heat conducting 
gas the conservation equations are assumed again to model the motion of a gas in one spatial 
dimension and time. Applying the conditions of expression (2.6) to equations (2.1) through (2.3) 
then yields 


5p 

dt 


, dp 
+ u^- + 
ox 



= o, 


(2.20) 



du 

+ p ‘% 


dp + d fyjdu ' | 
dx dx^ dx) 



du 
dx_ ’ 


(2.21) 



de 

+ pu Tx 


du d 
+ — 

dx dx 



(2.22) 


And assuming the Stokes’ relation holds 


X = 



(2.23) 


equations (2.21) and (2.22) become 
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(2.24) 


du , du _ dp 4 5 
P 3t pM 3^ “ _ 3x + 33x 


du 


de 


de 


du , 3 f , 3 


P ^ + pM ^ = -^ + ^J + 3 P U 


4 3 3m 


(2.25) 


Next, expressing the conservation equations (2.20), (2.24) and (2.25) in terms of the 
specific volume v = - rather than the density, and applying the subscript notation for partial 

p 


differentiation yields: 


v t - vu x + uv x = 0 , 

4 

u l + uu x +vp x --v([iu x ) x = 0, 

4 2 

e t + ue x + vpu x -v(k T x ) x - -p v (u x ) = 0. 

Now, from the equations of state (2.10) and (2.11) 

c„ 


e = pv ■ 


R 


Recall from thermodynamics that 
which yields 


R = c p~ c v > 


(2.26) 

(2.27) 

(2.28) 


e = 


pv 

Y - 1 ’ 


(2.29) 


where 



Combining equations (2.28) and (2.29) then produces the following form of the energy equation: 

■ [ipv) t + «C P v )x\ + v P u x ~ v i kT x) x ~ \^(u x ) 2 = 0 . (2.30) 

Equations (2.26), (2.27) and (2.30) are used to study shock jump conditions in an unsteady, one- 
dimensional, viscous, heat conducting, perfect gas. 
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3.0 NONSTANDARD ANALYSIS AND HYPERREAL NUMBERS 


3.1 Construction of a Hyperreal Number System 

A construction of a nonstandard or hyperreal number system *R is outlined in this section. 
The construction is roughly analogous to the construction of the real numbers R using 
equivalence classes of Cauchy sequences of rational numbers. For the hyperreal numbers, 
however, the equivalence relation is based on the rates of convergence of sequences of real 
numbers as well as the limits of convergence. The present notation and development essentially 
repeats that of Hurd and Loeb [24], Excellent elementary introductions to nonstandard analysis 
may be found in Henle and Kleinberg [25], Keisler [20,26] and Hoskins [27], while more 
advanced developments may be found in Lindstrom [28], Stroyan and Luxemburg [29], Davis 
[30], and Arkeryd, et al. [31]. 

The goal of a construction of a hyperreal number system is to produce a linearly ordered 

field (*R, +, ', <) which contains an isomorphic copy of the linearly ordered field of real numbers, 
but also contains infinitely large and infinitely small elements. Let Abe the set of natural numbers 
1, 2, 3, .... And let R s be the set of all sequences of real numbers indexed by N. Every element 
r e R s many be written as: 

r = ( r \> r 2’ r 3’ •••} = (3- 1 ) 

Natural componentwise operations of addition and multiplication may be defined for all r, s e R s 
via 

r®s = {(/„)} © {($„)} = {(r n + s n )}, (3.2) 

and 

r ® s = {(>■„)} ® {(*„)} = {(r n -s n )} = {(r n s n )}. (3.3) 

With these binary operations, the resulting 3-tuple is a commutative ring with the identity: 

1 = ( 1 , 1 , 1 , 1 ,...}, 

and the zero 

0 = ( 0 , 0 , 0 , 0 ,...}. 
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As an example of the algebraic operations defined by equations (3.2) and (3.3), consider the 
sequences 

r = { 1, 3, 5, 7, ... } and 5 = (2, 4, 6, 8, ... } . 

Then r ® s and r 0 s yield: 

r®s = (3, 7, 11, 15, ...}, 
and 

r®s = (2, 12,30,56, ...}. 

All algebraic concepts used in the subsequent development may be found in Hungerford [30], 

The space of all real sequences with the above addition and multiplication does not form a 
field, because of the existence of non-zero zero divisors, for example: 

(1,0, 1,0, 1,...} 0(0, 1,0, 1,0,...} = {0,0, 0,0,...} = 0. (3.4) 

This problem is repaired by introducing an equivalence relation ~ on R s and defining algebraic 

and order structures on a quotient space R s /~ to yield a linearly ordered field (*R, +, ', <). The 
equivalence relation is defined using the concept of an ultrafilter. 

Definition 1. Let / be a nonempty set. A filter on / is a nonempty collection U of subsets of / 
having the following properties: 

a. The empty set is not contained in U: 0 g U, 

b. If A, B e U , then A n B e U, 

c. If A e U and IziBziA then B e U. 

Definition 2. A filter is called an ultrafilter if for any subset A of I either A e U or its 
complement A c = I- A e U. 

Definition 3. If/ is an infinite set, the collection of subsets 7 7 / = {A<z.I\I — Ais finite} is a filter 
called the cofinite filter on I. More generally, set A is cofinite with respect to / if / - A is finite. An 
ultrafilter U on / is free if it contains Fj . 
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A free ultrafilter U cannot contain any finite set, or else the free ultrafilter U must contain the 
empty set and violates Definition la. The existence of free ultrafilters is not obvious and not 
trivial. Their existence follows from a form of the axiom of choice and is highly nonconstructive. 
Henle and Kleinberg [25] show how to obtain the elements of an ultrafilter on N using transfinite 
induction. The following theorem guaranteeing the existence of free ultrafilters is proved in Hurd 
and Loeb [24, p221]: 

Theorem 4. Free ultrafilters exist on any infinite set I. 

Now that it is known that free ultrafilters may be shown to exist on any infinite set, select 
any free ultrafilter TJ on the set of natural numbers N. U will contain all the cofinite sets of N. A 
relation — t/ is defined on the set of all real sequences of numbers, R s using the ultrafilter U. To 
simplify the notation the relation is written as ‘ ~ .’ 

Definition 5. If r, s <e R s , then r ~ s if and only if {n e N\r n = s n ) e U. 

Lemma 6. The relation ~ is an equivalence relation on R s . 

The proof of Lemma 6 is given in Hurd and Loeb [24, p4]. Since the relation ~ forms an 
equivalence relation on R s , R s is cut up into disjoint subsets or equivalence classes. Each 
equivalence class contains all the sequences of real numbers that agree on infinite indexing 
subsets of N which are contained in the free ultrafilter U. Following Definition 5, if 

r = (100,-5,7, 16, 16, 16, ...} and 5 = (-18, 11, 14, 15, 16, 16, 16, ...}, 
we have r ~ s , because the set {n e N\r n = s n } = (5, 6, 7, 8, ... } e U is cofinite with respect 
to N. The equivalence relation ~ is now used to define the nonstandard or hyperreal numbers. 

Definition 7. Let *R denote the set of all equivalence classes of R s induced by ~ : 

*R = R/~. 
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The equivalence class containing a given sequence s e R s is denoted by [.v], Therefore, if r ~ s in 


R s then 


[r] = [ 5 ] in *R. 

The elements of *R are called the nonstandard or hyperreal numbers. 


From the construction of *R via the free ultrafdter U it may be shown that either 
{1,0, 1,0, 1, . . . } or {0, 1,0, 1,0, . . . } but not both, will be identified with the zero element [0] 
in *R. Therefore, the ultrafilter removes the problem of non-zero zero divisors in the 
commutative ring based on the set of all real sequences, R s . Since the selection of the free 

ultrafilter U determines which sequence { 1, 0, 1, 0, 1, ... } or {0, 1, 0, 1, 0, ... } is identified 
with zero, the construction of the nonstandard numbers depends on U and does not appear to be 
unique. Uniqueness of the hyperreals turns out to be a deep and subtle issue. It has been proved 
that if the continuum hypothesis holds, the construction of the hyperreal numbers *R is 
independent of the selection of the free ultrafilter in the following sense: all such quotient fields 
constructed from a free ultrafilter based on N are isomorphic ordered fields. This means that one 
realization of *R is algebraically equivalent to any other realization of *R. However, the 
continuum hypothesis is not decidable within the Zermelo-Fraenkel axioms of set theory. If the 
continuum hypothesis is taken as an axiom, the hyperreals are unique; on the other hand without 
the continuum hypothesis, the uniqueness of the hyperreals is not decidable, Goldblatt [33]. 

Next to produce a linearly ordered field (*R, +, ', <), the algebraic and order relations must 
be defined. 


Definition 8. Let [r], [5] e *R. Then define (using equations 3.1 - 3.3): 

a. [r] + [s] = [r@s] = [{(r„ + $„)}] , 

b. [r]-[s] = [r][j] = [r®s] = [{{r n s n )}], 

c. [r] < [5] if and only if {n e N \ r n < s n } e U , and 

d. [r] < [5] ifandonlyif [r] < [5] or [r] = [5], 
where > and > are defined in a similar way. 
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The construction of *R together with the operations of Definition 8 yields a linearly ordered field: 


Theorem 9. (*R, +, <) is a linearly ordered field. 

The proof of Theorem 9 is simple; it is outlined in Hurd and Loeb [24, pp5-6]. 

Notice that the real numbers R from a linearly ordered subfield of (*R, +, ', <). This may be seen 
by introducing a mapping 

* : R — » *R, (3.5) 


to embed R into *R. 


Definition 10. If re R define the mapping of expression 3.5 by 

*(r) = *r=[{r, r, r, ...}] e *R. 

Theorem 11. The * mapping of Definition 10 is an order preserving isomorphism of R into *R. 
The proof of Theorem 11 is a straight forward verification of the relevant definitions. 

Definition 12. If A c: R, then A* is the set of all elements *a, where a e A . R* is the set of 
standard numbers in *R. 

Definition 13. Denote a hyperreal number by s (rather than [,s] ). 

a. A number s e *R is infinite if n < ,v for all standard natural numbers n; 

b. A number s e *R is finite if |.v| < n for some standard natural number n; 

c. A number s e *R is infinitesimal if .v| < - for all standard natural numbers n. 

At this stage, true infinitesimals may be identified in *R which do not exist in R. For 

example, 
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s = {1, 1/2, 1/3, 1/4, ...} = 



is an infinitesimal since given any non-zero and finite hyperreal number r = { ( r ) } , which is also 
a real number, s < r no matter how small r is. This is true, because if n 0 = [ 1/r] , the smallest 
natural number which is greater than or equal to 1 /r, then 

j« e N | - < r | = {n 0 , n 0 + 1, n 0 + 2, ... } e U . 

Note that the infinite sequence defining s is {(1 /«)} . Therefore, *R is much larger than R, 
containing infinitely small numbers which can be manipulated like finite numbers. Moreover, 
consider the hyperreal number 

n = {1,2, 3, 4,...} = {(«)}• 

This number is larger than any finite number p = {(p, p, p, and is therefore an infinite 
hyperreal number because the set {n e N \ n > p) e U . Note that 

Q+l = {2, 3, 4, 5, 6,...} = {(«+!)}, 



2 2 2 2 2 2 

Also note that Q“ = {1 ,2 ,3 ,4 ,...} = {(n ) } > Q gives another distinct infinite number. 
Notice that each finite hyperreal number is either equal to a real number or else it is 
infinitesimally close to a real number. 

Definition 14. A finite real number which is infinitesimally close to a finite hyperreal number is 
called the standard part of [a] , written st (a) or just a. 

For example if a = {( 1 + 1 / n)} , then st (a) = {(1,1,1,...)} = 1. 
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Definition 15. Let x,ye *R . 

a. x and y are infinitesimally close if x-y is an infinitesimal, written x « y . The monad 
of x is the set m (x) = {y e *R| x ~y} . 

b. x andy are finitely close if x-y is an finite. 

As implied by the examples above, a hyperreal number system contains many sets which 
are not subsets of the real numbers. The subsets of a hyperreal number system *R may be 
classified into two types of subsets: the internal and the external subsets. The internal subsets of 
are the elements of standard subsets of R embedded into *R by the * mapping of Definition 10. 

Definition 16. A subset A* c= *R is called internal if there exists a sequence { (A n ) } of subsets 
of R such that 

x = { (x n ) } e A* if and only if { n e N \ x n e A n } e U. 

A set 3 c; *R is external if it is not internal. 


The external objects contained in *R are the objects which do not result from objects in 
standard analysis, while the internal objects are precisely the objects which result from extending 
standard analysis via the * mapping. 


Definition 17. Let {(f n )} be a sequence of functions defined on a set A a R : 


fn- 


Then an internal function F 


is defined by the assignment 


F : A*->*R, 

x= {(x 75 )} e A* —> F(x) e *R, 


where A* ci *R and 


F(x) = {(/„)}•{(*„)} = {{f n (x n ))}. 
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Example 18. As an example of an internal function, consider the sequence of functions 


{/«(*)} = i I " 


n i + «v. 


An internal function F(x ) results for each hyperreal number {(x )} : 


F(x) 


\ n i + «x/ 


2 2 


If the sequences of real functions used to obtain an internal function F(x) are 
componentwise differentiable and integrable, the internal function will be differentiable and 
integrable via 

c (36) 

and 



( x - V 

J A(x)c/x = < 

a 

\fA x ) d ' 1 

Va„ J 


(3.7) 


where x = {(x n )} and a = {(<?„)}. 

Moreover, calculus over the hyperreal numbers is straight forward and is analogous to 
calculus over the real numbers. As examples, the derivative at x = a, a e R, of a differentiable 
function fix) is defined as 


dx 


= st 


r */(x)~*/(a) 
v x - a 


where x ~ a and x ^ a ; and the Reimann integral is given by the hyperreal sum 


Q 


J Ax) dx = 


= st 


£ *M k ) *x k 


a Vjt = 1 J 

where Q = [(«)] . See Keisler [20,26] and Henle and Kleinberg [25] for a development of the 
basic calculus in terms of the hyperreal numbers. 
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The * mapping has an important property that is key to applications of nonstandard 
analysis to standard real analysis. Suppose a logical sentence, denoted by O , is defined to be a 
formula in a first-order language for the real numbers. The truth of a sentence ® is conserved 
under the * mapping in a strong sense: *® is true if and only if ® is true. In practice this result 
means that theorems in real analysis are proved if nonstandard versions of the theorems are 
proved; that is, standard results follow from nonstandard proofs. This conservation property of the 
* mapping is called the transfer principle and follows from a theorem in mathematical logic on 
ultrafilter extensions due to Los. A proof of the transfer principle is developed in Hurd and Loeb 
[24], 


3.2 Nonstandard Generalized Functions 

In the early development of nonstandard analysis, hyperreal techniques were applied to 
the classical generalized functions, Robinson [15], with application to quantum theory, Kelemen 
and Robinson [34], Robinson noticed that any regularizing sequence of smooth functions used to 
represent the Dirac delta measure could be replaced by a single function with hyperreal 
parameters. Intuitively a single function represents the limit of the regularizing sequence. 
Consider the following delta sequence, 


S n ( X ) 


n 1 +nx 2 ' 


(3.8) 


Then letting the limiting parameter n be any infinite hyperreal number n Q , a smooth 
nonstandard function acting like the Dirac delta results: 


S(x) = - • 
n 


Q 

2 2 ' 

i + n V 


(3.9) 


This nonstandard delta measure behaves exactly like its classical counterpart: taking the standard 
part of the integral yields a functional producing the pointwise evaluation, 


00 00 

lim f s t1 (x)<\>(x)dx = st [-• 

J n 


D. 


y— oo 


n i + qV 


(| ){x)dx 


<K0), 


(3.10) 
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for an arbitrary test function (j) e D(R) . More generally, the right-hand-side of equation (3.10) 
may be used to define nonstandard linear functionals on spaces of internal test functions 
(|> e £>(*R). 

A nonstandard Dirac delta measure may be defined for any standard function g which 
satisfies 


00 

J g(x)dx = 1 , 

—00 

by considering the internal function 

S(x) = Cl ■ g(Q • x ) , 

where Cl is any infinite hyperreal number, Laugwitz [35]. 


Example 19. To see how the internal function of equation (3.9) acts like the Dirac delta, perform 
the integration: 


00 

J:- 


Q 


n i + fiV 


(j )(x)dx 


<K0) 


8 ( 0 ) 


n 


- 8 ( 0 ) 


Q 

2 2 

1 + Cl x 


dx 


(f)(0) 


n 


arctan(Qx)! 


where s = s(Q) is any infinitesimal such that Qs is an infinite hyperreal satisfying 


8 ( 0 ) 

- 8 ( 0 )’ 


Hence, 


arctan(Qs) 


n 


00 



—00 


— • () )(x)dx 

1 + Clx 


m 

n 



Then taking the standard part of the integral yields 

st(<|>(0)) = (f)(0) , 

which shows that the internal function of (3.9) acts like the delta measure. 

The power of the linear theory of generalized functions comes from the fact that many 
distinct regularizing sequences of functions may be applied to produce unique functionals, such as 
the delta measure. It is well known that the linear theory of generalized functions may be 
extended to a nonstandard theory of generalized functions, which reproduces the basic results of 
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the standard theory, Laugwitz [21] and Richter [36]. In the standard theory the delta measure is a 
unique functional resulting from the abstract limits of regularizing sequences; on the other hand in 
the hypcrrcal setting, the functional whose standard part produces the delta measure can be 
represented by uncountably many (internal) hyperreal functions with distinct microstructures. 
Such internal functions are called pre-distributions, Palmgren [37]. The significance of the 
nonstandard delta measures is in their utility; differentiable representations of such generalized 
functions may be written down explicitly and manipulated like standard differentiable functions. 
The focus of the present work is to exploit smooth representations of nonstandard generalized 
functions to derive shock wave jump conditions. No attempt is made in the present report to 
outline a nonstandard development of the linear theory of distributions. 

The analysis of shock wave jump conditions requires the use of hyperreal Heaviside 
functions. Like the delta measure, the Heaviside function can be represented by uncountably 
many differentiable internal functions with distinct microstructures. To introduce a nonstandard 
Heaviside function, let Sj and s 2 be two infinitesimal numbers where at least one of them is not 

zero, and where Sj « s 2 and Sj < s 2 . An internal Heaviside function H(x) is defined as follows: 


H(x ) = 


0 if v < a, 
h(x ) if Sj <x < s 9 , 

1 if *>s 2 


(3.11) 


where 

0 = *0 = ( 0 , 0 , 0 , 0 ,...}, 
h(x) = *h(*x) = {(h n (x n ))}, 
and 

1 = *1 = ( 1 , 1 , 1 , 1 ,...}. 

(Here and in the following development the *-notation is suppressed.) Moreover, h(x) is assumed 
to be a differentiable function, not necessarily monotonically increasing, which satisfies 

h{ gj) = 0 and h{ e 2 ) = 1. (3.12) 

As an example, the following function can be used for h(x) when x is in the open infinitesimal 
interval (s 1? s 2 ) . 
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Example 20. Figure 1 shows the standard Heaviside function 


^std( x ) 


Jo if x < 0 
{ 1 if x > 0 


and the nonstandard Heaviside H(x ) function given by combining equations (3.13) and (3.1 1): 

0 if x < 0 

~ ~ cos — if 0<x<s, (3.14) 

1 if x > s 

where s is a positive infinitesimal. All the plots of the hyperreal functions are normalized on the 
finite interval [0, 1] by setting x — > x/s . Notice that the standard Heaviside function has a jump 
at x = 0 but the nonstandard Heaviside function is continuous [0, s] . Also notice that H(x ) is a 
differentiable function on (0, s) since h(x) is differentiable. The derivative of II(x) is the 
internal function 

0 if x < 0 

sin — it 0 < x < s , 

_ 8 _ 

0 if x > s 

which is a pre-distribution for the Dirac delta. For all real test functions cj) e D(R) it may be 
shown that 





following a calculation similar to that performed in Example 19. An important property of H'(x ) 
is 
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J H'(x)dx = ^ • J sin[7rQx](/x = 1, 


e(Q) 


where Q is any infinite number and s = s(f2) is any infinitesimal number such that Qs = 1 . 
Figure 2 shows the derivative of the nonstandard Heaviside function H'(x ) of equation (3.14) on 
the normalized interval [0, 1] . 

Another important property is that the support of H'(x) is the closed infinitesimal interval 
[0, s] . Suppose \|/(x) is standard function which is discontinuous at x = 0 with right and left 
limits \|/(0+) and \j/(0-) respectively. Then because s > 0 , the integral satisfies 


st 


J H'(x ) • \\f(x)dx 


00 


V(0+). 


This means that for a function such as \\ /(x) with discontinuity at x = 0 , the derivative of II(x) 
picks up \|/(0+). By selecting s 1 and s 2 in equation (3.11), another nonstandard Heaviside 
function G(x) may be found such that G'(x) will have the following property 


st 


J G'(x) ■ \\f(x)dx 


00 


a\)/(0-) + (1 - a)\|/(0+), 


where 0 < a < 1 . 

8 £ 

In particular, let s be a positive infinitesimal and Sf = s 2 = - in equations (3.11) 
and (3.12) to define a differentiable Heaviside function: 


K(x) = 


0 if *<-§ 


1 1 

- - - cos 
2 2 


n(s/2 +x) ~ 
8 


• j, 8 8 

it - “< X < - , 

2 2 


1 


if 


Then K'(x) has the following property 
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St 


( \ 

00 

J K'(x) ■ \\f(x)dx 

V— GO ) 


= |(Y(0-) + Y(0+)). 


3.3 Multiplication of Nonstandard Generalized Functions 

The Schwartz theorem showing that generalized functions cannot be multiplied is a 
manifestation of a basic structural limitation of algebraic mathematics. Following Rosinger [5], 
consider the mathematical structure (3, °, D) consisting of a collection 3 of functions defined 
on R which may exhibit a finite number of localized discontinuities. The standard Heaviside 
function (with the STD notation suppressed) 


H(x) = 


0 

1 


if x < 0 
if x > 0 ’ 


(3.15) 


is contained in 3 . Let ° and D denote the standard multiplication and differentiation operators 

acting on the set 3 . 

From equation (3.15), notice that 

77™ (x) = ( 77°77°77°...°77)(x ) = 77(x), (3.16) 

for all m > 1 . Then applying the derivative operator in the sense of generalized functions to 
equation (3.16) yields 

mH' n ~\x)DH(x ) = DH(x), 

which implies that 

m H m ! (x) = 1 , 
and 

H(x) = (3.17) 

m 

for all m> 1 . But equation (3.17) contradicts the definition of the standard Heaviside function 
(3.15). Therefore, the mathematical structure (3, °,D) is over-specified and discontinuity, 
multiplication and differentiation cannot be combined in the algebraic steps of a single derivation. 
The over-specification of (3, °, D) may be repaired by requiring that 
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H ( x) = (H°H)(x) * H(x) . 

Recall that a hyperreal number system is a field; hence, hyperreal numbers may be multiplied. 
This implies that the hyperreal functions may be multiplied on intervals where they are defined. 
Note that for any Heaviside function defined by equation (3.11), 

0 if x < s j 


H\x) 


h (x) if Sj < x < s 7 


(3.18) 


1 


if x>s- 


which implies that 


iz (x) ^ H{x ) , 


2 

because h (x)*h(x) for s 1 <x<s 7 . That is, giving the Heaviside function H(x ) a fixed 

microstructure removes the algebraic over-specification which leads to the contradiction of 
equation (3.17). Also note, that by taking derivative of both sides of equation (3.18) 


H(x)H'(x) = \ 


0 if x < 8 j 

h{x)h'{x) if s 1 <x<s 9 . 
0 if x > s 2 


For example, equations (3.11) and (3.13) yield 


0 if x < 0 


H(x)H'(x ) = 


n . 
— sin 

4s 


7TX 

8 


1 - COS 


7TX 

8 


if 0 < x < s 


0 if x > s 


Therefore, H(x)H'(x ) has a clearly defined hyperreal range assuming s is a given nonzero 
infinitesimal. Equation (3.11) generalizes for all differentiable nonstandard Heaviside functions to 
the result that 

H m (x)*H(x), 

for all m > 1 , removing the contradictions obtained by setting 

//”(x) = H(x) . 
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Example 21. Although the product of two smooth pre-distributions exists and is a smooth 
nonstandard function, such a product will not yield a nonstandard distribution. To see this, recall 
the internal function of equation (3.9) representing the delta measure and the consider the 
functional acting on any standard test function cj> e D( R) ; 


00 

J S(x) • 8(x) • () )(x)dx 

—00 



Q 


“|2 


1 + Q 2 x 2 J 


• c) )(x)dx , 


(3.19) 


where Q is any infinite hyperreal and s = s(fi) is any infinitesimal hyperreal such that Qs is an 
infinite hyperreal. Performing the integration yields 

e(Q) 


1 Q I 2 


J[iT 


2 2 

+ fix 


• cj ){x)dx 


(f)(0) 

f 

r Q 1 

2 

n 

J 

2 2 

L 1 +Q x“ J 


dx , 


-s(fi) 


where 


1 


1 Q l 2 


n 1+Q 2 x 2j 


() ){x)dx 


' 777X1 + 2 k ' arctan < n *)) 

n z 1 + Q x 


e(fi) 


-s(Q) 


Then evaluating the limits of the integral 


— f • q2 (j T7 + 777 ’ arctan(Qx) 

n V2 1+Qx IL1 


s(Q) 


-s(fi) 


<t>(0) o 2 
n 


1 2 

1 + Q 8 


+ 


£)• 


Notice that the first term of the right-hand-side of the integral satisfies 


m 


Q. 


n 


Qc_) _ c|)(0) Q 


2 2 

1 + Q 8 


2 1 

” ^i + Q£ 


m 1 

2 8 


n 


while the second term may be written 

*9 ■*■(.£) 


d>(Q) 

2n 


■ Cl . 


Therefore, the functional for the square of delta measure with the pre-distribution of (3.9) 
produces 


J 8(x) • 8(x) • <\>(x)dx = J - 


1 Q 


n 1+qV 


•<K x)dx * <|>(0) • f — 2 ■ - + (3.20) 

v n £ zn 
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which is an infinite hyperreal number. And since the integral of the square of the Dirac delta is an 
infinite number, the standard part of the integral is infinite. Hence, the square of the delta measure 
does not yield a generalized function; and, the specific hyperreal number produced for 8-8 
depends on the choice of the pre-distribution used to represent the delta measure. 


Example 22. The result of equation (3.20) and its interpretation for the square of the delta 
measure changes dramatically if the functional of equation (3.19) is assumed to act on a space on 
nonstandard internal test functions D(* R) instead of a space of standard test functions D( R) . To 
see this, choose any internal test function (j) e D{* R) such that 

(f)(0) = s 2 , 

where Q is any infinite hyperreal and s = s(Q) is any infinitesimal hyperreal such that Qs is 

2 

infinite and Os” is finite. Integrating equation (3.19) then yields 


GO 

J S(x) • S(x) • cf )(x)dx « (f)(0) ■ ( — 2 ■ ^ ^ • o] « ~ 2 + ^ • n • s 2 , 

Ti 


(3.21) 


and taking the standard part of the integral produces the finite value: 

f ~ \ 


st 


J S(x) • 8(x) • (f )(x)dx 


= ^st ( n. e 2 ). 


V-00 J 

Hence, the square of the delta measure (3.21) yields a nonstandard object which is analogous to a 
Colombeau nonlinear generalized function. The specific finite number produced for 8 • 8 
depends on the choice of the pre-distribution used to represent the delta measure as well as the 
choice of the nonstandard test function. 


Example 23. As a second example of a product of generalized functions, consider 8 • H . The 
multiplication of the delta measure with the Heaviside function occurs in the analysis and 
computation of shock waves if the governing equations of motion are written in nonconservative 
form; 8 • H comes from the nonlinear terms, such as the convective acceleration term uu x , for a 

discontinuity in the field variable, u . In the standard theory of generalized functions, the product 
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8 • H does not yield a well defined functional, because the results depend on the representations 
of 8 and H . To see this behavior, choose the following pre-distributions for the delta measure and 
the Heaviside function: 



Q 

exp(Qx) + exp(-Qx) ’ 


H(x) = 


0 if x < 0 

exp (-fir) if 0 < r < s , 

1 if x > s 


where Q is any infinite hyperreal and s = s(Q) is any infinitesimal hyperreal such that Qs is an 
infinite hyperreal. Now, consider the functional acting on any standard test function (j) e D( R) : 


00 

J S(x) • H(x ) • (j )(x)dx 

—00 


00 



0 


n ^~. exp( ~ n ^ ' • Ux)dx 
exp(Qx) + exp(-f2x) 


Integrating equation (3.22) yields 


(3.22) 



s(Q) 

Q( 1 - exp(-Qx)) a r Q( 1 - exp(-Qx)) 

exp(Qx) + exp(-Qx) n n \ J exp(Qx) + exp(-Qx) 

o 


where 


00 



0 


f2( 1 - exp(-Qx)) 
exp(Qx) + exp(-Qx) n } ' 


4 >( 0 ) 


n 


{ 2arctan((exp(Qx)) + ln(l + exp(2Qx)) - 2Qx)} 


s(n> 

o 




Hence, 


st 


J 8(x) • H(x) ■ §(x)dx 


= 0.27936 • c()(0) 


GO 


(3.23) 


If other pre-distributions are selected for 8 and H , the numerical value multiplying c|)(0) in 
equation (3.21) will, in general, be different. In the numerical approximation of shock problems, 
the use of nonconservative forms of the equations of motion produce shock wave speeds that 
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depend on the numerical scheme and not the physics of shock propagation. Nonconservative 
numerical simulations of one-dimensional shock problems may be found in Baty [38], Such 
numerical approximations are analogous to the result of equation (3.21); the specific numerical 
values are fixed by the choice of the pre-distributions for 8 and H . 

Example 24. Again, consider the product 8 • H , but now assume 8 = H' . Then for any pre- 
distribution representing the delta measure, H' defined on the interval (0, s), choose hyperreal 
numbers Q, s = s(Q) so that Q is infinite, s is infinitesimal and Qs is an infinite hyperreal. 
Then for any standard test function cj> e D( R) , the functional equation for 8 ■ II implies: 

00 00 

J 8(x) • H(x ) • (j )(x)dx = J H'(x ) • H(x) ■ §(x)dx, 

—00 —00 

and, 

00 oo s(Q) 

fS (x) ■ H(x) ■ <\>(x)dx = f — [\ ■ H (x)l • (j)(x)r/x ~ cj)(0) f —[\-H(x)]dx . (3.24) 

J J dx^2 J J dx v 2 J 


Therefore, 


st 


J S(x) • H{x) ■ (j )(x)dx 


= 2 ^( 0 ) • 


(3.25) 


J 

Equation (3.25) is consistent with the results obtained by the standard theory of generalized 
functions. This special case of the product 8 • H may be computed using only discontinuity and 
multiplication and is not over-specified, since 8 • H is not differentiated. The differentiation of 
the Heaviside function is avoided by applying integration by parts. To this end, consider the 
standard version of equation (3.24): 


00 00 CO 

f 8(x) • H{x) ■ cj ){x)dx = f — [\ ■ H (x)l • (j)(x)dx = f — [\ ■ H(x)\ §(x)dx , (3.26) 

J J dx v 2 J J dx v 2 J 

— 00 —00 —00 

2 

which assumes that H(x) = H (x) . Equation (3.24) then implies: 
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J • H ( x )) ty(x)dx = • ^^-§(x)dx = ^ ■ 4>(0) , 

-00 0 

producing the same result as equation (3.25). This example is analogous to integrating the 
governing equations of motion modeling shock wave propagation in conservative form: the 
predicted shock wave speed will be independent of the numerical approximation. 

Notice that the standard calculation performed in Example 24 to reproduce the result of 
equation (3.25) depends on the order of the algebraic and analytic steps. Changing the order of the 
operations can produce an over-specified problem with ambiguous results. This demonstrates the 
power of applying nonstandard analysis to problems involving products of general functions: the 
mathematical inconsistencies are removed and the problem may often be reduced to elementary 
calculus. 
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4.0 INVISCID SHOCK WAVE JUMP CONDITIONS 


4.1 Specific Volume, Velocity and Pressure Jump Conditions 

Consider the motion of a normal shock wave propagating through a quiescent, perfect gas. 
On either side of the shock wave the specific volume v , velocity u , pressure p , and entropy 5 , of 
the flow are assumed to be constant. The goal of this section is to derive the jump conditions for 
specific volume, velocity and pressure across a normal shock assuming that the governing 
equations of motion are expressed in nonconservative form. The jump conditions for specific 
volume and pressure are then applied in the following section to derive the jump condition for 
entropy. Since the values of the field variables are constant on either side of the shock wave, the 
shock wave does not accelerate, so that the shock speed is constant. Hence, the characteristics for 
the shock wave are straight lines in the space-time domain. Along the characteristic lines the 
governing equations reduce to total or ordinary differential equations, so the field variables v , u , 
and p , may be assumed to have the following form of equation (1.2) across the shock wave: 


v(5) = v ( +[«), (4.1) 

k(S) = Ut + MKtt), (4.2) 

P($) = Pi +[p]m), (4.3) 

where £, = x + ct is a characteristic line, c is the shock speed, and where [())] = <j) ; . - <\>[ , with 


the subscripts r and / referring to the right (downstream) and left (upstream) conditions across the 
shock, respectively. The functions //(£, ) , K(t ) , and L(c) are assumed to be smooth (internal) 
nonstandard Heaviside functions as described in Chapter 3. Each one of these Heaviside functions 
is taken to have its jump contained on the same infinitesimal interval (0, s) . 

Now substituting equations (4.1) and (4.2) into the continuity equation (2.17) yields: 

c[v]H' - ( V/ + [v]H)[u]K' + (uj + [u\K)[v\H' = 0 , (4.4) 


where X' 


d.V 


Defining u = u + c and collecting terms in equation (4.4) produces 


(U[+ [u]K)[v\H' - (v i + [v]H)[u]K' = 0, 


which implies: 
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M H' = [u]K' 
V/+ [v]/7 Uj + [u]K' 


Integrating both sides of equation of (4.5) then gives: 


v, 

H = - M + alM + K 


[v] y[u] 


(4.5) 


where a depends on the constant of integration. Applying the fact that H,K = 0 for £, = 0 , 


yields: 


and hence, 


V/[»] 

MV 


H = 


vj[u] 

Mu, 


(4.6) 


Recalling that H, K = I for c = a , then implies: 

V ,M = 

Mu, 

And combining equations (4.6) and (4.7) yields, 

H = K, 

so that H and K are the same function defined on the infinitesimal interval (0, s) , i.e., H and K 
have the same micro structure. Notice that equation (4.7) is equivalent to the classical shock wave 
jump condition for conservation of mass: 

P ,u, = P r u r . 

Next, substituting equations (4.1), (4.2) and (4.3) into the momentum equation (2.18) 

yields: 

c[u]K' + (Uf+ [u]K)[u]K' + (vj+ [v]77)[p]Z' = 0. (4.9) 


(4.7) 

(4.8) 


Applying H = K from equation (4.8) produces the separable equation: 


(ii/ + [i/]A)[w]if' 
v , + [v]A 


+ [p]L' = 0. 


(4.10) 


Then combining equations (4.7) and (4.10) implies: 
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(4.11) 


Integrating equation (4.11) gives: 


j^K' + [p]L' = 0. 


L = -r4^+P, 


[v][p] 

where p is a constant of integration. Applying the property that L, K 
P = 0 ; while applying L, K = 1 for £, = s , implies: 


0 for £, 


[uV 


- 1 . 


[v]|>] 

Equation (4.13) is equivalent to the Prandtl relation: 

U r U/ = |>][v] = 


(4.12) 
0, gives 

(4.13) 


Pr-Pl 
P r ~ Pi 


Whence, equations (4.12) and (4.13) show that the microstructure for the jump functions K and L 
is the same: 

K = L. (4.14) 

Finally, by combining equations (4.8) and (4.14) it follows that the microstructure for the 
Heaviside functions for the specific volume, velocity, and pressure jump conditions are coincident 
across an inviscid shock wave: 

H = K = L. (4.15) 

An important observation should be drawn from the analysis producing equation (4.15). 
The algebraic and analytic operations shown in equations (4.4) through (4.13) appear to follow 
closely well known calculations such as those performed by Colombeau and Le Roux, [10] and 
Salas and Iollo [18] yielding the same result. In these works, however, the objects being 
manipulated are nonlinear generalized functions (functional equations) written in a symbolic 
notation similar to that used in advanced calculus; these objects are not smooth functions and do 
not support all of the operations of ordinary algebra and calculus; the multiplication of singular 
generalized functions is accomplished via a weak equality called association. In contrast to such 
calculations, the objects manipulated in equations (4.4) to (4.13) (and indeed, in the following 
sections of this report) are smooth nonstandard functions. Therefore, the algebraic and analytic 
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manipulations are the standard operations for smooth functions extended to smooth internal 
functions. 


Example 25. Notice that the derivation of equation (4.15) by using the governing equations in 
nonconservative form requires manipulating products of the pre-distributions, //(£,) , K(t ) , and 
/.(£,) and their derivatives. Since these functions are differentiable the mathematical structures 
applied in this analysis are not over-specified and will produce the same results obtained by 
analyzing the problem in conservative form. For example, an expression for the shock speed c 
may be derived using equation (4.7): 

y /M = | 

[v]«/ 


This relation implies 

V/|>] 

u, = — — 

[v] 

and since u l = u l + c , it follows that 


c = ^(v^j-^v]). (4.16) 

Equation (4.16) is the shock speed obtained from the equations of motion in conservative form 
(which is also the physically observed shock speed). 


4.2 Entropy Jump Condition 

To derive the entropy jump condition, as in section 4. 1 assume a solution for the entropy s 
in the form of equation (1.2) across the shock wave 

s(S) = Si+[s]S&), (4.17) 

where £, = x + ct is a characteristic line, c is the shock speed, and where [ 5 ] = s r - s l , with the 

subscripts r and / referring to the right and left conditions across the shock, respectively. 
Combining equations (4.2), (4.3), (4.14), and (4.17) with the energy equation (2.15) yields: 

(Pi + b]E0(c[5]y + (u, + MOONS') = 

c v (c[p]H' + (u ! + [u]H)[p]H' + yip, + [p]H)[u]H') 
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which implies 


( Pj + [p]H)(ui + [u]H)[s]S ' = c v ((u,+ [u]H)[p]H' + y(p l + [p]H)[u]H') 
Equation (4. 18) then produces: 


MS' = cJ l 


Recalling that equation (4.7) gives 


Pi + [/;]/// '-it i + \u\hJ 


V/[«] 


u l ~ r i ’ 

[v] 

implies that equation (4. 19) may be written: 

\pW 


= c. 


+ y 


Pj + [p]W \\j + [v]H' 


[yW 


(4.18) 


(4.19) 


(4.20) 


Integrating equation (4.20) then produces: 

[s]S = c v log((/? z + [p]H) • ( v z + [v]//) Y ) + C , 
where C, is a constant of integration. Now, applying the property that H = 0 for ^ 

C = -c v log(p z v}); 


(4.21) 
0 , gives 

(4.22) 


while applying H = 1 for £, = s , yields, 

M = <M«g((l + ^ • («3) 

Equations (4.21), (4.22), and (4.23) then yield the jump function for entropy across an inviscid 
shock wave: 


[s] 


[p] 


F^log i + • 1 + —77 


Pi 


[v], 


(4.24) 


Notice that equation (4.24) may be derived from the thermodynamic equation for the 
entropy, equation (2.19): 

5 = s z +[s]S = c v \og((p,+ \p\H)-(v l +[v]H) 1 ). 


Hence, 
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S = CylogCp/v/) + c v log((l + ^ H ) ■ (i + ’ ( 4 - 25 ) 

where 


s i = c v log O/v/) . 

Equation (4.25) shows that the jump function for entropy is a nonlinear function of H . If 
the Heaviside function H is assumed to be smooth and strictly monotonically increasing on the 
infinitesimal interval (0, s) , then the entropy jump function exhibits a peak in the shock layer: s 
increases from s l at 0 to a peak and then decreases to s r at s such that s l < s r . To see this, 
differentiate equation (4.25) to produce: 

d^ = ([pVPiW Y (M/V/)#' 

Cv ’i +m/p,w Cv ' i+([v]/v,)/r '■ > 

To find an extreme point, set 


cE 

d| 


= 0. 


(4.27) 


Since H is a smooth, strictly monotonically increasing function, equations (4.26) and (4.27) give 

[ %+([v]/v / )//) + Yv J (l +([p]/pi)H) = 0, 

Pi v / 


which implies 


[£] +y m +(r+ 1 ) [£]m ff( y , o. 


Pi 


Pi v / 


(4.28) 


Recall from gas dynamics that pressure and specific volume jump terms and may be related 
through the Rankine-Hugoniot equation for a perfect gas: 


[v] = ~2[pVPi 
V/ 2y + (y+l )[p]/pi’ 


(4.29) 


Thompson [23]. Now, combining equations (4.28) and (4.29) produces: 


[p] / -2[p]/pi \ 

Pi V2y + (y+l )[p]/p/ 


+ (y+l) 


\pV ~2[p\/Pi 
Pi ^2y + (y + 1 )[p]/p 


H(U = 0 . 


which reduces to 
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and hence, 


2r + (y+l)^-2y-2(y+l)M//(5 ) = 0, 

Pi Pi 

m p ) = \ (4.30) 

And from the assumptions on H , a point ^ e (0, s) exists statisfying equation (4.30). This 
derivation yields the following result: 

Theorem 26. If H is a smooth, strictly monotonically increasing, internal Heaviside function, 
the entropy jump function of equation (4.25) admits a maximum value in the shock layer (0, s) . 

From both theoretical and experimental work, it is known that H is a strictly 
monotonically increasing function, Muntz and Harnett [39]. Furthermore, the existence of the 
entropy peak in the shock layer is well known and has been derived using a viscous flow 
assumption, Morduchow and Libby [17] and using Colombeau’s nonlinear generalized functions, 
Salas and Iollo [18]. The present analysis demonstrates that the entropy peak is caused by 
combining a nonlinear expression with a jump condition and does not depend on the viscosity of 
the fluid. 

To plot examples of the normalized entropy jump in the shock layer, the smooth internal 
Heaviside function is assumed for the specific volume: 

H(f) = ^ ~ • arctan , (4.31) 

where Q is an infinite hyperreal and s is an infinitesimal hyperreal such that Qs is an infinite 
hyperreal. On the interval [0, s] the Heaviside function of equation (4.31) satisfies: 

H( 0) « 0 and H( s) * 1 . 

Figure 3 shows equation (4.31) on the normalized jump interval of [0, 1] with Q = 100 . Now, if 
the pressure and specific volume jump terms [ ] , of equation (4.24) are written as functions of the 
Mach number, M : 
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(4.32) 


and 


[p1 = 

Pi Y + 1 


1 ) 


[v] 

v l 



(4.33) 


with the upstream flow variables set to p l = 1 and v l = 1 ; the entropy jump function may be 
plotted. Figure 4 shows the variation of Sit,/ s) on the normalized jump interval [0, 1] for a 

7 

perfect gas flow with y = - for the normal shock wave Mach numbers M = 1.5, 2.0, 2.5, 3.0 . 

These entropy jump functions shown in Figure 4 exhibit peak values as required by 
Theorem 26. Notice that the peaks all occur at t,/e = ^ . This is the result of the choice of the 
Heaviside function in equation (4.31). As shown in the proof of Theorem 4.26, the peak of the 
entropy occurs where the jump function obtains the value H(^ p ) = ^ . Since, 

H( 1 / 2 ) = l ~, 

the peak occurs at the mid point in the jump interval. Recall that the inviscid analysis does not 
determine the specific Heaviside function for the field variables; only the relationship between the 
field variable jump functions is found. Hence, different representations of the jump function for 
specific volume will yield different entropy shapes and entropy peak locations. The key point of 
the theorem is that a peak must occur for any strictly monotonically increasing jump function 
assumed for the specific volume. 


Example 27. An analysis similar to that used to derive the entropy jump condition shows that 
across the shock wave: 


pis t + us x ) ^ 0 . 

Suppose that equation (4.34) is zero across a shock wave. Recall from Chapter 2, that 

pi s t +us x) = c v(p t +u Px + yp u x)’ 


so that 


(4.34) 
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p(s t + us x ) = 0, 


implies: 

Pt+up x + yP u x = °- ( 4 -35) 

Substituting equations (4.2) and (4.3) into equation (4.35) and applying equation (4.17) produces, 

((«/+ [u]H)[p] +y(Pi+ [p]H)[u])H' = 0, 

which may be written: 

(«/l>] + JPi[u] + (Y+ \)[u][p]H)H' = 0. (4.36) 

Equation (4.36) implies that either 

H' = 0, (4.37) 


or 


Ui[p]+YPi[u]+(y+ l)[u][p]H = 0, 


(4.38) 


must hold. Equation (4.37) cannot be hold because H is not constant; hence, equation (4.38) must 
be considered. By applying the values on the ends of the jump interval, £, = 0 and £, = s , to H , 
equation (4.38) yields two expressions: 

«/l>] + YP/I>] = °> (439) 

and 

u,[p] + yp 1 [u] + (y+ l)[u][p] = 0. (4.40) 

Then combining equation (4.39) with equation (4.40) produces: 

[u][p] = 0. (4.41) 

Equation (4.41) implies that either the velocity or the pressure is constant across the shock wave, 
which is physically incorrect (a basic contradiction). The upshot is that none of the inviscid shock 
wave jump conditions may be derived in an unambiguous way by assuming 

s t + us x = 0 . 
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5.0 VISCOUS SHOCK WAVE JUMP CONDITIONS 


5.1 Jump Conditions for a Viscous Non-Heat Conducting Gas 

Consider the motion of a normal shock wave propagating through a quiescent, viscous, 
perfect gas. For the analysis of this chapter, the notation defined in section 4.1 will be applied. 
The goal of this section is to derive the jump conditions for specific volume, velocity, pressure, 
and entropy across a normal shock assuming that the governing equations of motion are expressed 
in nonconservative form. For the present analysis the coefficients of visocity p and thermal 
conductivity k are assumed to be constant. To derive these jump conditions, assume solutions for 
the field variables v, u, p , and s in the form of equation (1.2) across the shock wave: 


V(5) = v,+ [v]tf(S), 

(5.1) 

«(4) = «/ + [«]£(£), 

(5.2) 

p(£) = Pi+[p]L(Z>), 

(5.3) 

s($) = SJ+WS&). 

(5.4) 


The functions //(c) , ATc) , L(t) , and S(t) are smooth nonstandard Heaviside functions with 
their jumps contained on the infinitesimal interval (0, s) . 

Applying equations (5.1) and (5.2) with the continuity equation (2.26), a relationship 
between the Heaviside functions for specific volume and velocity may be determined. This 
analysis is identical to the inviscid derivation shown in equations (4.4) through (4.8) which 
produces: 

H = K. (5.5) 

Hence, the specific volume and the velocity have the same microstructure across the shock wave. 

Next to develop a relationship between the Heaviside functions for the specific volume 
and the pressure, equations (5.1), (5.2) and (5.3) are substituted in the momentum equation (2.27): 

c[u~\K' + («;+ [u]K)[u]K' + (V;+ [v]F/)[p]2/ fe|p(v / + [v]H)[u]K" = 0. (5.6) 

Applying equation (5.5) to equation (5.6) produces: 

c[u]H' + (uj+ [u]H)[u]H' + (vj + [v]F/)[p]T r - ^p(v ; + [v]H)[u]H" = 0, 
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which implies, 


(u, + c)[u] + [wl H 4 

v, + [v]n "■W-jtoVr-o. 


(5.7) 


Combining u t = u t + c with equations (4.7), and (5.7) yields: 

2 

l^H' + [p]L' -^[uW = 0, 
or 


TTrr [A] Tit \p~\ T I — A 


(5.8) 


where (3 = -p . Integrating equation (5.8) with respect to the characteristic variable gives: 


H' 


r»l JT = I p]_r +r 
PM P [«] 


(5.9) 


here ^ is a constant of integration. Applying the boundary conditions H,H',K= 0 for q = 0 , 
shows that the constant of integration is zero: 

C = 0. (5.10) 


And applying the boundary conditions H, K = 1 and H' = 0 for £, = s , then implies: 

[uf = -MM, (5.11) 

which the same condition found using the inviscid momentum equation, (4.13). Equations (5.9) to 
(5.11) give: 


H' + 


M 

PM 


H = 


M J 
PM ’ 


(5.12) 


which is a linear, non-homogeneous ordinary differential equation. Integrating (5.12) with respect 
to the characteristic variable £, and applying the boundary condition 77(0) = 0 then produces: 


% 

H(£) = 0jexp(0{q-^})Z(q)r/q , 
o 


where ® = 


M 

PM 


and hence establishing the relationship between H and L . 


(5.13) 
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Continuing, another equation relating the Heaviside functions H and L may be derived 
using the energy equation (2.30) for a viscous, heat conducting, perfect gas. Recall that: 

pv = RT. (5.14) 

Notice that assuming the field variables satisfy the functional form ( 1 .2) it follows that: 

(pv), = c(pv) x . (5.15) 

Then combining equations (5.14) and (5.15) with equation (2.30) implies 

1 2 

— [ u (pv) x + vpu x -r\v(j>v) xx -MUxY = (5.16) 

k 

where r\ = - . From equations (5.1), (5.2), (5.5), and (4.7) it follows that 
K 

Ur 

u = -v. (5.17) 

v l 

From (5.17) equation (5.16) then becomes 

a(pv) x + pu x -p(pv) xx -$(u x ) 2 = 0, 

and hence, 

a(p x v+pv x )+pu x -i~[(p xx v + 2p x v x +pv xx )-p(u x ) 2 = 0. (5.18) 

, 1 ~ u i 

where a = • — . 

y - 1 v i 

Next, applying equations (5.1), (5.2), and (5.3) with H = K and to equation (5.18) then 

yields 

a [p]L'(vj + [v]H) + a [v]H'(p l + [ p]L ) + (p l + [p]L)[u]H' 

~ T l[p]C"(v / + [v]//) - 2r[(p)(v)L' H' - r\(pj + [p]T)[v]//" - p[u] 2 (//') 2 = 0. (5.19) 
To express equation (5.19) in terms of H , equation (5.12) implies: 

[p]L = ^[u]H' + [p]H , 

[p]L' = P[u]H" + [p]H', (5.20) 

[p]L'' = ^[uW' +\PW ■ 

Substituting equations (5.20) into (5.19) then yields: 

a(P[u]/7" + [p]H')(V/ + [v]H) + (a[v]TT + [u]H')(pi + P [u]H' + [p]H) 
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-q(p [u\H'" + [p]H")(v l + [v]H)-2r)[v]H'($[u]H" + [p]H ') 

-^(p l + V[u]H' + [p]H)[v]H"-p>[u]\H') 2 = 0. (5.21) 

Combining like orders of derivatives in equation (5.21) produces: 

-r|P 

+ (ap[w](v z + [v]//)- r) [p](v l + [v]77) - 2qp [u][v]H' - rj (p, + p [u]H' + [p]H)[v])H" 

+ a[p]H'(v/ + [ y]H) + (a[v] + [u])H'(p l + p [u]H' + [p]H ) - (2q[p][v] + p = 0 , 
so that 

-riP [u\H"'(vi + [v]H) 

+ (aP[w]v / ^ri[p]v / ijri[v]/> / )if" + (ap[i/][v] -2r\ \p]\v])HH" - 3qp [v]\u]H'H" 

+ i^[p]v l + a[v]p,+ [u]pj)H' + (afi[u][v]-2r][p][v])(H') 2 

+ (a b][v] + M^)(// 2 )' = 0, 

and collecting the terms for ( HH , 

-r|p [u]H"'(vi+[v]H) 

+ (cl^[u]v ! - r\[p]v r\[v]p ^H" + (ap[w][v] - 2r\[p][v])(HH')' - 3^^[v][u]H' H" 

+ (a[p]v, + a[v]^ ; + [«]p ; )F + (a[p][r] + *)(^)' = 0. (5.22) 

Now for the case of a non-heat conducting gas, set k = 0 which implies q = 0 . 
Equation (5.22) then reduces to 

a P \u]V/H" + ap [u][v](HH')' + (a [p]v z + a [v]pj + [u^p^H' 

+ (a [p][v] + [ -^)(H 2 )' = 0. (5.23) 

Recalling equation (4.7) and 

1 ~ u i 

a = • 

Y - 1 v, 

imply that equation (5.23) may be written: 
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(5.24) 


U$MH" + («/[p] + JPi[u])H' + p [u]\HHy + y -^[p][u](H 2 y = o. 

Integrating the last expression yields: 

U[$[u]H' + (u^p] + y p l [u\)H + p [u] 2 HH' + - ^ [p][»]/7~ = k, (5.25) 

where k is a constant of integration. The boundary conditions //, //' = 0 at c, = 0 then imply 

k = 0, (5.26) 

while the boundary conditions H = 1 and H' = 0 at £, = s give: 

M/Lp] +Y^/[m] !>][«] = 0- (5-27) 


Example 28. Equation (5.27) may be shown to be true by rewriting this coefficient using equation 
(4.7) as: 


t£] +Y M + Lj_L[£lM = o 

Pi vj 2 Pl v t 

Combining equation (5.28) with the Rankine-Hugoniot equation (4.29) then yields 


l+r(: 


-2 


+ 


Y+Uhif 


-2 


which implies 


as claimed. 


2y + (Y+1 )\pVp/ 2 p, V2y + (y+ \)[pVp 


2y + (Y+ 1)Lp]/^/-2y-(Y+ \)WPi = 0, 


= 0 , 


(5.28) 


Applying equation (5.26) to (5.25), the ordinary differential equation for the Heaviside 
function becomes: 


and 


(w/PM + fi[u] H)H' + {ii,[p] +yp,[u] + ^-^[p][u]H}H = 0, 


Pm(i + ^h) h ' + + = °- 


(5.29) 


Finally, recasting equation (5.29) as an integral equation yields: 
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1 +^H 


Mr 


Pi 


Pi v l 




(5.30) 


where C, is a constant of integration. 

To compute the Heaviside function for the specific volume, equation (5.30) was solved 
numerically. Since the left-hand-side of (5.30) is zero for \ = 0, the constant of integration 


becomes zero, C, = 0 . The coefficient 


M. 

PM 


appearing on the right-hand-side of equation (5.30) 


contains the viscosity of the gas and is a very large number which implies that the jump will occur 
over a very small spatial interval. Once equation (5.30) was solved for H{f ) , equation (5.20) was 
applied to solve for Z(^) . 

Figures 5 and 6 compare the computed Heaviside functions for specific volume and 
pressure for the Mach numbers of 1.5 and 3.0. All calculations performed shown here assume 



These Figures show that the microstructure associated with the specific volume and the 


pressure is distinct and is a function of the Mach number of the flow field. A key result of the 
viscous analysis is that the equations of motion contain enough physical information to specify 
the jump functions H, K and L . Figure 7 shows the specific volume jump function in terms of the 
Mach numbers 1.5, 2.0, 2.5 and 3.0; while, Figure 8 shows the pressure jump function for the 
same Mach numbers. 


5.2 Jump Conditions for a Viscous Heat Conducting Gas 

The jump conditions for a viscous, heat conducting, perfect gas follow from the analysis 
of section 5.1 with the minor modification that the coefficient of thermal conductivity k is 
assumed to be a non-zero constant. Recalling that the functions //(c) , K{t ) , and, L(c) are 
smooth nonstandard Heaviside functions with their jumps contained on the infinitesimal interval 
(0, s) for the specific volume, velocity and, pressure, respectively. The jump conditions for these 
field variables are governed by equations (5.5), (5.13), and (5.22): 

H = K, (5.31) 
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= 0jexp(0{q-^})Z(q)^, 


(5.32) 


and 


-^[uWivt+MH) 

+ (ap[n]v / -ii[p]v / -ii[v] j p / )//" + (a^[u][v]-2^[p][v])(HH'y -3^[v][u]H'H" 

+ (ab]v z + a[v]^ ; + [«]^)F + (a[p][v] + *)(^)' = 0, (5.33) 


where 

^ \p] o 4 k 4 1 ~ u i 

O = -jp- , p = -]x , r| = - , and a = • - . 

P[n] 3 R y-1 v l 

Equation (5.33) may be integrated once to yield a second order ODE for the specific volume jump 
function: 

7\ 2 


-ilP[n]i/"(v / + [v]/E)-iiP[n][v](//') 

”( r l[vlP/ - ap[«]v / + r||j?]v / - ap[n][v]//+ 2r|[p][v] 77)77' 


+ 


([u\pi + a[v]p/ + a\p]vi)H+^([u] + 2a[v])H 2 = k. 


(5.34) 


By applying the boundary conditions 77, 77' = 0 at c = 0 , 

k = 0, 

while the boundary conditions 77 — 1 and 77' = 0 at £, = s are satisfied since equation (5.34) 
reduces to equation (5.27). The boundary conditions for equation (5.34) then become: 

77( 0) = 0 and 77( s) = 1 . (5.35) 

Once equation (5.34) is integrated to determine the specific volume Heaviside function, 77, the 
pressure Heaviside function is computed via equation (5.20): 

Ls m H ' +H , 

[p] 

To compute 77 using equation (5.34), the coefficients are grouped so that the physical 
parameters associated with the pressure and specific volume jumps across the shock wave may be 
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expressed in terms of the Mach number by applying equations (4.32) and (4.33). The visocity and 
heat conduction parameters are introduced in the Prandlt number. The Prandlt number enters into 
equation (5.34) via the ratio of p and p : 


g_p_£ _3_ _ 7 _3_ y(y - 1 ) _ 3 


y 


(5.36) 


p R 4\i R 4p y(y-l) 4 (y-1 )Pr' 

Equation (5.30) for a viscous non-heating gas, shows that the thickness of the shock depends on 
[nl 4 

the parameter O = ' containing the viscosity [1 = - p . To remove this parameter from the 
P [ u ] 3 


governing equation, the following change of variables is introduced: 

q = -®^, 

so that 

_d_ r y = _l_d_ 
dq K> Od§‘ 


(5.37) 


(5.38) 


For derivation in this section the Prandlt number is fixed to be Pr = - . So equation (5.36) 


implies: 


8 = 


y-1 

Then applying equations (5.36), (5.37), and (5.38) to (5.34), the governing equation reduces to: 

2 


8—H" + M hh" + (( — + s') ^ + S^ET + 8^tzJ(7/' 
Pi Pi v, vvy -1 y Pi V Pi v i 


) 


+ 


- — — + sl — —HH' + 
2(y- 1) y Pi v. 


y-i y Pi v 


H- 


1 y + 1 [p][v] 

y-i 2 Pi v/ 


(H) =0. (5.39) 


To find the specific volume Heaviside function, //(q) , for a viscous, heat conducting gas, 
equation (5.39) must be solved numerically. To perform such calculations, the boundary 
conditions of (5.35) were replaced with initial conditions. The initial data may be specified from 
the physics of the flow field requiring that the jump function satisfy an isentropic compression as 
the flow enters the shock wave. In practice, this amounts to fixing the initial data as: 

7/(0) = 0 and 77(0 ) = 3(M, 9), (5.40) 

where & is a very small value. Numerically it was observed that as /V7 1 or if M was fixed and 

& — > 0 the computed jump function II(t) becomes symmetrical about an inflection point. 
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Figures 9, 10, 11 and 12 show Heaviside functions for specific volume, pressure, entropy 
and temperature, computed for a viscous, heating-conducting gas as a function of the Mach 

7 

numbers 1.5, 2.0, 2.5 and 3.0. As in the above section, all calculations shown assume y = - . The 

temperature Heaviside function was computed from the specific volume and pressure jump 
functions using the equation of state (2.11). Figures 13 and 14 compare the microstructure of the 
field variables for the Mach numbers 1.5 and 3.0. The entropy for a viscous, heating-conducting 
gas exhibits a peak as was predicted via the inviscid analysis. The microstructure peaks observed 
in the entropy and the temperature are the result of jump functions combined with nonlinear 
equations of motion. 


5.3 Comparison of Shock Wave Jump Conditions with Experimental Data 

To check the shock wave microstructure predicted by equation (5.39), the resulting 
specific volume jump function, //(£, ) was compared with the density measured across a rarefied 
shock wave by Muntz and Harnett [39]. In the Muntz-Harnett experiment, a supersonic helium 
flow with a Mach number of 1.59 was used to generate a normal shock wave. The resulting 
normal shock wave microstructure was approximately 2.2 cm wide, a distance which represented 
-3-4 mean-free-paths for the HE molecules. The data produced in the experiment was also shown 
in reference [39] to be in good agreement with a computational solution of the Navier- Stokes 
equations. 

For the present comparison, the Prandlt number and ratio of specific heats were specified 
for helium as: 



Moreover, density was nondimensionalized using the following functional form: 


with 


p* = 


P-P/ 
Pr - Pi 


(5.41) 


Pi = 1.92 x 10 8 g/cm 3 and p r = 3.51 x 10 8 g/cm 3 . 

Substituting the computed specific volume Heaviside function into equation (5.41) then yielded: 
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l/( Vl +[v]H)- Pl 

p* = . 

P r ~ P / 

Figure 15 compares of the measured shock wave microstructure with that predicted by equation 
(5.39). Overall the predicted density (specific volume) jump function is in good agreement with 

the measured density. The main variation observed occurred near p* = 1 . The shape of the 
density jump function depends on the initial data, (5.40). As the initial data approached zero the 
predicted results produced were in good agreement with the experimental data. The predicted 

jump function matched the measured jump condition closely around p* = 0 . 

The micro structure computed for the other field variables for normal helium shock wave 
are very similar to the jump functions shown in Figures 10, 11 and 12 of Section 5.3 for the 

3 7 

M = 1.5 case with Pr = - and y = - and, as such, are not shown. 

4 ' 5 
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6.0 SUMMARY AND CONCLUSIONS 


This report derived shock wave jump conditions for a one-dimensional compressible gas 
using nonstandard analysis. Differentiable internal Heaviside functions were introduced and used 
to model the microstructure of the flow field across a normal shock wave governed by the 
classical conservation equations. For flow variables with their jumps defined on the same 
infinitesimal interval, it was shown that the equations of motion written in nonconservative form 
produced unambiguous relations between the various nonstandard Heaviside functions. The 
analysis specified the microstructure across the shock wave for the cases in which the equations 
of motion contained a physical model of the shock region. For an inviscid shock wave, the 
analysis only produced the relations between the field variable Heaviside functions, because the 
equations of motion contained no physical information about the shock region. On the other hand 
for a viscous shock wave, the analysis predicted the micro structure across the shock wave. The 
entropy peak observed across a shock wave was predicted by both the inviscid and the viscous 
analyses and is the result of combining jump conditions with nonlinear field equations. 

The present report was motivated by the authors’ interest in understanding how to add 
physical information to specify products of generalized functions contained in differential 
equations modeling physical problems with well defined regions of solution discontinuity. The 
main conclusions of the report are as follows: 

1. Nonstandard analysis may be used to identify differentiable internal Heaviside 
functions defined on infinitesimal intervals that reproduce the standard jump 
function results from distribution theory. 

2. The derivative of a differentiable internal Heaviside function yields a pre- 
distribution of the Dirac delta measure. 

3. Examples of differentiable internal Heaviside functions and pre-distributions 
of the delta measure were applied to compute nonstandard versions of products 
of generalized functions. The products of such generalized functions were 
shown to depend on the specific pre-distributions used to represent the 

2 

product. Higher order singular products such as 5 in general do not yield 
standard generalized functions for functionals defined on spaces of standard 
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test functions; although standard generalized function-like objects may result 
for functionals defined on spaces of nonstandard test functions. 

4. Application of differentiable pre-distributions of the Heaviside function and 
the delta measure in products of generalized functions removes the algebraic 
over-specification and the associated mathematical ambiguities. 

5. While mathematical ambiguities are removed by applying differentiable pre- 
distributions, additional physical information must be given to find problem 
specific pre-distributions for a given product of generalized functions. 
Products of generalized functions may be defined consistently via nonstandard 
analysis; however, physically meaningful products of generalized functions 
must be specified from the physics of the problem under investigation. 

6. The relations between the internal Heaviside functions for a one-dimensional 
compressible inviscid gas were derived unambiguously from equations of 
motion in nonconservative from by using differentiable pre-distributions to 
represent the products of generalized functions. 

7. The problem specific internal Heaviside functions for a one-dimensional 
compressible viscous gas were found (computed) unambiguously from 
equations of motion in nonconservative form by using differentiable pre- 
distributions to represent the products of generalized functions. 

8 The shock wave jump conditions were computed directly by solving for the 
specific microstructure of the Heaviside functions of the flow parameters 
instead of the flow parameters themselves. The flow microstructure was 
defined on the characteristic lines as a traveling wave to transform the partial 
differential equations of motion to ordinary differential equations. The 
application of nonstandard analysis with the characteristic geometry of the 
flow produced a significant simplification of the problem, since the resulting 
ordinary differential equations are easier to solve than the equations of motion 
used in classical analyses. 
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8.0 FIGURES 



Figure 1: Comparison of standard and nonstandard Heaviside functions. The normalized 
jump interval is the unit interval [0,1], 
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Figure 2: The derivative of the nonstandard Heaviside function shown in Figure 1 
plotted on the normalized jump interval. 
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Figure 3: A smooth Heaviside function based on the inverse tan function plotted on the 
normalized jump interval with Q =100. 
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Figure 4: The inviscid entropy Heaviside function as a function of the Mach number 
plotted on the normalized jump interval. 
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Figure 5: The viscous specific volume and pressure Heaviside functions for the Mach 
number 1.5, plotted on the normalized jump interval. 
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Figure 6: The viscous specific volume and pressure Heaviside functions for the Mach 
number 3.0, plotted on the normalized jump interval. 
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Figure 7: The viscous specific volume Heaviside function as a function of Mach number 
plotted on the normalized jump interval. 
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Figure 8: The viscous pressure Heaviside function as a function Mach number plotted on 
the normalized jump interval. 
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Figure 9: The viscous, heat conducting, specific volume Heaviside function as a function 
of Mach number plotted on the normalized jump interval. 
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Figure 10: The viscous, heat conducting, pressure Heaviside function as a function of 
Mach number plotted on the normalized jump interval. 
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Figure 11: The viscous, heat conducting, entropy Heaviside function as a function of 
Mach number plotted on the normalized jump interval. 
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Figure 12: The viscous, heat conducting, temperature Heaviside function as a function of 
Mach number plotted on the normalized jump interval. 
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Figure 13: The viscous, heat conducting, field variable Heaviside functions for the Mach 
number 1.5, plotted on the normalized jump interval. 
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Figure 14: The viscous, heat conducting, field variable Heaviside functions for the Mach 
number 3.0, plotted on the normalized jump interval. 
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Figure 15: The comparison of the density micro structure computed with Equation (5.39) 
with the density micro structure measured for a Mach number 1.59, normal, rarefied 
helium shock wave. 
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